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ABSTRACT. The weak gravitational lensing of high redshift type la super- 
novae has the potential of probing the structure of matter on galaxy halo 
scales. This is complementary to the weak lensing of galaxies which probes 
structure of larger scales. There are already several organized searches for 
these supernovae being carried out for the purposes of cosmological parameter 
estimation. A method is proposed for extracting from future supernovae data 
information on lensing and the structures responsible. This method utilizes 
the correlations between SN luminosities and foreground galaxies. Simula- 
tions of the lensing and uncertainties will be presented. It is found that with 
a hundred supernovae or more at z ^ 1 or larger, significant measurements 
of the mass, shape and extent of dark matter halos could be made if they 
contain a significant proportion of the matter in the universe. 



1 Introduction 

Two of the most active and 
promising fields of observational cos- 
mology today are the search for weak 
gravitational lensing by large scale 
structure (LSS) and the search for 
high redshift supernovae (SNe) for 
the purposes of determining the lu- 
minosity distance-redshift relation. 
These two subjects are discussed 
elsewhere in this proceedings (Hook 
1998, Schneider 1998 and Bartel- 
mann 1998, these proceedings). Here 
I will discuss another use for high 
redshift SNe, and one that is comple- 
mentary to weak lensing measures of 
LSS. 

Others have addressed the lens- 
ing of SNe in the past. These in- 
vestigations have concentrated on ei- 
ther microlensing (Linder, Schneider 
& Wagoner 1988) or the effect lens- 
ing has on uncertainties in cosmolog- 
ical parameter estimation (Frieman 
1997, Wambsganss et al. 1997, Kan- 
towski 1998, Holz 1998). Here I ad- 
dress the problem of detecting lens- 
ing and recovering information about 
the lensing structures from future SN 
data. 

A little review of the SN searches 
will be necessary. The most impor- 
tant point here is that after correc- 
tions type la SNe have been shown 
empirically to be good standard 
candles with peak luminosities of 
Mg ~ -19. A variance of AMy t B — 



0.12 mag in the corrected peak mag- 
nitudes has been achieved using V 
and B bands ( Reiss, Press & Kirsh- 
ner 1996 ) and a AM B ^ 0.17 mag 
using one color (Hamuy et al. 1996). 
It should be noted that these num- 
bers are only approximations of the 
true intrinsic variation in peak lumi- 
nosities because they are based on 
the statistics of just a few SNe (20 
and 18, respectively). The technique 
for correcting the peak luminosities 
may improve in the next few years as 
more SNe are found and additional 
information such as spectral features 
are incorporated. At present over 100 
type la SN have been observed, most 
of them below z ~ 0.5. Searches can 
now reliably discover about 10 SNe 
in a night of observing (~ 3 deg 2 ). 
For more details on this subject see 
Hook in these proceedings and refer- 
ences there in. 

In the next section some motiva- 
tions for considering the lensing SNe 
will be given. In section 3 simulations 
of the lensing are described and in 
section 4 these simulations are used 
in Monte Carlo calculations to esti- 
mate the accuracy of future experi- 
ments. 



2 Why Supernovae? 

It is reasonable to ask the ques- 
tion: Why try to detect the lensing of 
SNe when you can use galaxies which 
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are much more numerous and can be 
observed at higher rcdshifts? The an- 
swer is that the lensing of SNe and 
galaxies probe different size scales in 
the lensing structure. 

Unfortunately we do not know the 
intrinsic properties such as elliptic- 
ities and luminosities of galaxies at 
2^1. Not only do we not know 
the properties of individual galax- 
ies, but because of galaxy evolu- 
tion and the fact that galaxies are 
lensed we also do not know the in- 
trinsic statistical distribution of any 
observable properties that are af- 
fected by lensing. For example we do 
not know the rms ellipticity of un- 
lcnsed galaxies at z 2. As a re- 
sult galaxies taken individually can 
give us no information about lens- 
ing. To search for lensing we must 
look for correlations amongst the im- 
ages of galaxies - background and 
foreground or background and back- 
ground - which are assumed to not 
exist in the real galaxies.^] The ne- 
cessity of using correlations between 
galaxies makes studies of the lens- 
ing of galaxies insensitive to scales 
below roughly the mean separation 
between galaxies even though these 
scales probably dominate the lens- 
ing of individual galaxies. Noise in- 
creases this minimum scale to ,C 1'. 

Supcrnovae type la on the other 
hand are presumably the same ob- 
jects today as they where at any 
redshift excepting the possibility of 
evolution in progenitor properties. 
By observing low redshift SNe we 
can determine the unlensed statis- 
tical properties of these sources - 
i.e. their peak luminosity distribu- 
tion - and then take each high red- 
shift SNe to be an independent mea- 
sure of the magnification at a point 
which in general will be a larger sig- 
nal than the magnification averaged 
over some area. 

There is another reason why the 
lensing of SNe should give a higher 
signal. Matter on small scales is ex- 
pected to be centrally concentrated 
in clumps. The shear, -f(r), is a "po- 
lar" (a vector like object that rotates 
according to 28 instead of 8) so it has 
a direction. The influence of different 
mass clumps on a particular galaxy 



image can partially cancel. The shear 
will then add like a random walk. 
This is important for galaxy-galaxy 
lensing which attempts to measure 
the shape of galaxy halos. This tech- 
nique has been used to measure the 
mass within ~ 30/i _1 kpc, but does 
not put strong constraints on the 
shape of halos (see Hudson et al. 
1998 and references there in). In the 
weak lensing limit the magnification 
/i(r) — 1 oc S(r), where S(r) is 
the projected surface density of the 
clump. Magnification has no cancel- 
lations. 

Because the variance in the point- 
wise magnification is larger then the 
variance in the smoothed magnifi- 
cation it might be hoped that SNe 
will make up for their small numbers. 
This is the question I wish to address 
in the rest of this paper. 



3 Simulations of Lensing 

To simulate the results of possi- 
ble future observations a model for 
the lensing of point sources must 
first be made. In the weak lensing 
limit the magnification can be writ- 
ten in terms of the convergence, k(8), 
fi(8) ^ l + 2/e(0). Then the change in 
a sources' magnitude is 

Am = 2.51og(/j(0)) 2.17k(0). (1) 

In this paper I will only be inter- 
ested in the mass correlated with ob- 
servable galaxies (for other possibili- 
ties see Metcalf 1998). It will be as- 
sumed that this mass is in discrete 
halos surrounding each galaxy. PI The 
convergence due to one halo is 

Khirx) = -=T? T = K h -K b (2) 

^c(Zh,ZSN) 

The second term, Ej, accounts for 
the compensating effect of the back- 
ground density. This term is often 
left out of this expression when con- 
sidering the lensing by clusters be- 
cause it is small and shear maps are 
insensitive to a uniform offset. In the 
present case it will be important. The 



1 An exception to this is the case of 
strong lensing were the distortion in 
the image is very large. 



2 The term "halos" is perhaps a lit- 
tle over restrictive because what is 
really measured is the correlation be- 
tween galaxies and mass. 
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Figure 1. This is the distribution of convergences, n, for SNe at z = 1 in 
Q = 0.3 universes, both open and Q a + Q\ = 1. The vertical lines on the 
left sides of the plots demark the "empty beam" solutions as explained in the 
text. 



convergence calculated without this 
term is and E c is the critical den- 
sity. 

The value of depends on what 
distance is used to propagate the 
light between halos. Here the usual 
Robertson- Walker distance will be 
used and will be set so that 
(k) = and so the mean luminosity- 
redshift distance will agree with 
Robertson- Walker. This can be done 
experimentally. If the lensing is suffi- 
ciently weak the total magnification 
of a SN will be the sum of the mag- 
nifications of each halo. 
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(3) 



The dark halo of each galaxy is 
modeled with the density profile 



p(r) 



K 2 



i-wGr 2 



+ 1 



(4) 



This is just an isothermal sphere 
sewn to a power law at a scale length 
r s . The central velocity is taken to 
be related to the luminosity of the 
galaxy through the Tully-Fisher rela- 
tion, V c = VfiL/Lff. A Schechter 
function fitted to local galaxies is 
used for the luminosity function. 
This model assumes that the halos 
were formed sometime before the SN 
went off and that there has not been 
substantial merging since that time. 
For 7 < 2 a cutoff radius is usually 
needed so that the halos do not over 
fill the universe. This additional pa- 



rameter can be expressed as the total 
mass in halos relative to the critical 
mass, n^. 

The clustering of halos is done by 
running several relatively low reso- 
lution N-body simulations with dif- 
ferent random initial density fields. 
These simulations are stopped at dif- 
ferent redshifts and a group finding 
algorithm is used to find the loca- 
tions of halos. Halos of the form of 
equation (^) are then put in, random 
lines of sight are propagated through 
the simulation and the magnification 
of each SN is calculated using equa- 
tion §. 

Figure 1 shows the results of some 
of these simulations. It is clear that 
the medians of the distributions are 
all less than K = 0. Most galaxies are 
demagnified. It is also evident that 
s< 1 for the vast majority of SNe at 
this rcdshift. This justifies assump- 
tion that the lensing is weak. 

The minimum possible value of K 
occurs when the line of sight does 
not intersect any halo. This might 
be called the "empty beam" approx- 
imation for the angular size distance 
although in this case not all of the 
mass is in the halos so the beam is 
not truly empty. This minimum k 
is represented by the vertical lines 
in figure It can be seen that in 
A models the minimum re is smaller 
and the distribution extends to lower 
k than in open models. In addition 
it is more likely in a A model that 
there will be a close-encounter with 
a halo. A result of these two things 
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is that the variance of P(re) is larger 
in a A model than in an open model, 
all other things being equal. The left 
panel in figure |l] shows that the clus- 
tering of the halos does not make 
a great deal of difference in their 
lcnsing properties in these models. 
This means that the important scales 
for the lensing of point sources are 
,$ IMpc, a range not accessible in 
shear maps except in the atypical 
cases of clusters. 

There is information in probabil- 
ity distribution function of re (the re- 
p.d.f) not only on the background 
cosmology (which can be gained in 
other, better ways) and the density 
profile of the halos, but also on the 
nature of dark matter. I have im- 
plicitly assumed in this paper that 
the dark matter can be treated as a 
smooth density distribution for the 
purposes of lensing. If the dark mat- 
ter is in planet or stellar mass objects 
the re-p.d.f. will change drastically. 
The peak of the distribution will be 
taller and narrower and very close to 
the minimum re. For more on this use 
of the re-p.d.f. see Metcalf (1998). 

4 Simulations of Observations 

The goal here is to predict how 
well future observations of SNe could 
constrain parameters of the halo 
model. Or conversely, how much and 
what kind of data will be needed 
to make a good measurement of the 
halo parameters. 

The method I propose here is to 
cross-correlate properties of observed 
foreground galaxies with the cor- 
rected peak luminosity of the SNe. 
Using photometric redshifts and a 
trial Tully-Fisher relation a trial re 
for each SNe can be calculated. Af- 
ter a number of SNe have been ob- 
served the halo model can be ad- 
justed to best fit the data. The model 
outlined in the previous section has 
the following set of free parameters: 
{Vf, r s ,7, /3, f2^}, if L j is taken to be 
fixed. Some of these would normally 
be held fixed as others are adjusted 
to fit the data. 

The major sources of uncertainty 
are in the photometric redshifts, the 
scatter in the Tully-Fisher relation 
and the scatter in the intrinsic peak 
luminosities of the SNe. Of these the 
latter by far dominates the uncer- 
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Figure 2. A two dimensional his- 
togram representing the results of a 
Monte Carlo simulation for Q = 
0.3, n h = 0.18, tt A = 0.7 with 7 = 1 
halos. The shading represents how 
many realizations of the observation 
fell within that region. The contours 
contain %63 and %95 of the realiza- 
tions. 

tainty. 

Monte Carlo calculation is done 
to assess how well observations are 
likely to constrain parameters. A set 
of N SNe is created with their as- 
sociated lensing halos. Many realiza- 
tions of random errors are then gen- 
erated to make simulated data sets. 
For each data set a set of best fit- 
ting parameters is found. The spread 
of these recovered parameters about 
the set that was used as input is a 
measure of an experiment's uncer- 
tainties so far as the input model is 
realistic. 

Figures |^ and |i| shows the results 
of two such simulations. Here the 
constraints on the combination of Vf 
and r B , the scale size of the halos, 
are shown if 300 SNe at 2 = 1 are 
observed. The variance in the SNe's 
peak magnitude is taken to be cr m = 
0.12. There is a bias in the recov- 
ered parameters, but it is small in 
comparison to the uncertainties. The 
main difference between these two 
cases is in - 7 is not strongly con- 
strained. 

Simulations with 100 SNe at z = 
0.9 and a m = 0.16 still show a signif- 
icant detection of lensing and a con- 
straint on log(r s ) of <r = 0.44. Al- 
lowing the scale length to be corre- 
lated with the galaxies' luminosity, 
like for example r 3 oc \/L, increases 
the dispersion in re, and thus the sig- 
nal, because more mass is put into 
rarer bright objects. 
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Figure 3. The same as figure |^ only 
with Sift = 0.29 and 7 = 3. 



There are a number of large field 
of view CCD cameras coming on line 
in the next few years. Combining 
searches for lensing through shear 
with searches for SNe using these 
cameras would be mutually benefi- 
cial since the lensing of SNe provides 
information on scales not accessible 
in shear maps. 

It should also be mentioned that 
it may also be possible to detect 
the lensing of SNe without resort- 
ing to correlations between matter 
and light. This subject is addressed 
in Metcalf (1998). 



5 Discussion 

How many SNe are there? If the 
comoving type la SN rate is con- 
stant it is expected that there are 
500 - 2000 SNe/dcg 2 /yr with z < 
1.5. The larger end of this range is for 
the low Q models. The SN needs to 
be caught within about a one week 
(times 1 + z) window to be useful. 
This gives ~ 15 - 80 SNe/deg 2 be- 
low z = 1.5 at any given time. In ad- 
dition since the star formation rate 
goes up with redshift it stands to rea- 
son that the SN rate goes up as well. 
This could increase the SNe rate by 
as much as a factor of 5 (see Ruiz- 
Lapuente P. & Canal R., (1998), Sa- 
dat et al. (1998) and others for dis- 
cussions of the la SN rate). The SNe 
are certainly out there. 

The detection of type la SNe at 
z = 0.97 has already been demon- 
strated although this SN was not 
spectroscopically identified as a type 
la (Garnavich, et al. 1998). The spec- 
troscopic identification of the SNe as 
type la is the limiting factor in push- 
ing to higher redshift right now. In 
addition to high redshift SNe, low 
redshift SNe will be needed to ac- 
curately determine the distribution 
of SN peak luminosities and hope- 
fully to improve the techniques used 
to correct the peak luminosity. 

With the simple model presented 
here it has been demonstrated that 
with 100 or more SNe at z > 0.9 
something meaningful can be said 
about the nature of dark matter ha- 
los. Higher redshift would help of 
course, but even with z = 0.9, <r m = 
0.16 and N S n ~ 100 lensing can be 
detected if most of the mass in the 
universe is within a few Mpc of a 
galaxy and is made of small particles. 
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